Prediction and formation mechanism of triaxial superdeformed nuclei for A ~ 80 
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The three dimensional Total Routhian Surface (TRS) calculations are carried out for 64 nuclei 
between 70 < A < 90 to find triaxial superdeformed nuclei. Total of 12 nuclei were predicted to have 
triaxial superdeformation in which the neutron rotational energy plays a key role and the neutron 
shell energy plays additional role in the formation of triaxial superdeformed nuclei. 
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I. INTRODUCTION 

Usually, the shape of a deformed nucleus is supposed 
as a ellipsoid with small hexadecapole deformation. In 
nilsson model, the frequency of harmonic oscillator is de- 
scribed as 



u k = UJ [1 
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where £2 is the quadrupole deformation parameter and 7 
is the triaxial deformation parameter. Let the half axis 
of x, y, z are a, b, c, respectively, then from au3\ = &w 2 = 
CCJ3 , we get the relations between a, 6, c and £2, 7: 



£2 = 
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Rf [62 cos 7 + V3e2 sin 7 + 3] " 
Rf [£2 cos 7 — V3s2 sin 7 + 3] " 
Rf[3-2e 2 COS7]- 1 , 



(4) 



here / = [-2e| cos(37) - 9s| + 27] 1/3 and R is the radius 
of a sphere whose volume is equal to the ellipsoid volume. 
For example, when a : b : c = 1 : 1 : 2(axial ellipsoid), the 
£2 and 7 will take the volume of 0.6 and 0°, respectively. 
Such nuclei, i.e., c : a ~ 2 : 1 and a ~ b, is called axial 
superdeformed nuclei. While a =/= b, the triaxial deforma- 
tion 7 will not be zero. For instance, if a : b : c — 3 : 4 : 6, 
then (£2,7) will be (0.577,30°). Such nuclei which has 
large 7 deformation and large quadrupole deformation is 
called triaxial superdeformed (TSD) nuclei. From Eq.© 
we can see that 7 is sensitive to the difference of a and b. 

Until now, many data for axial superdeformed nuclei 
are accumulated but the data for triaxial superdeformed 
nuclei are quite few. Only five nuclei, 163 Lu, 165 Lu, 
167 Lu, 171 Ta, 86 Zr[l-5] were identified experimently to 
have triaxial superdeformation. the triaxial behavior of 
163 Lu has been further confirmed experimentally by the 
discovery of the wobbling mode |]|. However, the tri- 
axiality of most these TSD nuclei can only be explained 



by the theoretical calculations. Among these five discov- 
ered TSD nuclei, four are located in A ~ 160 region and 
only one is in the A ~ 80 region. The prediction of tri- 
axial superdeformed nuclei in A ~ 160 region has been 
done [6]. From that work several other nuclei besides the 
four discovered, were predicted to have triaxial superde- 
formation for proton configuration [660] 1/2. It is also 
pointed out that the shapes may co-existences for other 
curtain configurations and a nucleus may have triaxial 
superdeformation in different q.p. configurations. 

The TSD nuclei in A ~ 80 , 86 Zr, was discovered in 
1998. It does not seems like that TSD 86 Zr is a accident 
appeared in the A ~ 80 region. Other TSD nuclei must 
exist in this region. In this paper, we attempt to give 
the prediction of the TSD nuclei near A ~ 80 by TRS 
(Total Routhian Surface) calculations and their forma- 
tion mechanism . Section 2 includes brief description of 
the three-dimensional TRS theory which is used to de- 
termine the nuclear deformation; Section 3 includes the 
prediction of TSD nuclei in A ~ 80 region; Section 4 in- 
cludes the discussion of the formation mechanism of the 
TSD nuclei; Section 5 gives a summary. 



II. A BRIEF DESCRIPTION OF THE TRS 
MODEL 

The hamiltonian of quasi-particles moving in a 
quadrupole deformed potential rotating around the x- 
axis with a frequency lo may be written as 

H" = JT B . P . (£ 2 , £4, 7) ~ AAT + A(P + P+) - uj J x , (5) 

where H s p denotes the deformed hamiltonian of sin- 
gle particle motion, the second term on the right hand 
side is the chemical potential, the third term is the pair- 
ing interaction and the last term stands for the Coriolis 
forces. The modified-harmonic-oscillator (MHO) poten- 
tial with the parameters k and fi for the mass region 
taken from Ref.[7] is employed in the present calculation. 
The pairing-gap parameter is determined empirically by 
A = 0.9A O C ., and A o c . is taken from experimental odd- 
even mass difference [8]. As an approximation, we did not 
take the deformation and rotation dependence of pairing 
into account. 
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The total routhian surface, namely the total energy in 
the rotating frame as a function of the £2, 7, and £4, of 
a (Z, N) nucleus for a fixed quasi-particle configuration 
c.f. may be calculated by 

E c - L (e 2 ,e il i]uj) = E\ A (e 2 ,e 4: ,i) + E COII (e 2 ,e 4 ,-f,uj = 0) 

iGc.f. 

Where E\a is the liquid-drop model energy [9], -E^orr is the 
quantal-effect correction to the energy, which includes 
both the shcll[10] and pairing corrections[ll]. The col- 
lective rotational energy E TO t may be calculated micro- 
scopically as the energy difference between the expec- 
tation values of with and without rotation, by us- 
ing the wave function for the quasi-particle vacuum con- 
figuration^]. The last term of Eq.(2) is the sum of 
quasi-particle energies belonging to the configuration c.f., 
which generates the deformation drive. All of the term 
in Eq.(2) depend on the {Z,N) numbers which are not 
written explicitly. The equilibrium deformations of nu- 
cleus may be calculated by minimizing the total routhian 
energy of Eq.(2) with respect to £2, £4, and 7. Here we 
put the hexadicapole deformation £4 as a free parameter 
in order to get better results. 

In the real process of minimizing the total routhian, 
we minimize the total routhian as a function of £4 for 
each point (£2,7) and get two surface E cA -(e2,"f) and 
£4(^2,7) first. Then from the surface of E c - L (e2,j) we 
can find the minimum -E" mm an< ^ the corresponding £2 m in 
and 7 min . From the £4^2,7), ^2min and 7 min , we may 
determine the value of £4 m i n - From these steps above, 
we may find the equilibrium deformation £2 mm, £4min 
and 7 m in which possibly exists in the nucleus. 



III. THE PREDICTION OF THE TRIAXIAL 
SUPERDEFORMED NUCLEI 

Before the calculation of A ~ 80 nuclei, the TRS 
method was checked and compared with the discov- 
ered TSD nuclei, 86 Zr. Our calculated result, (£2,7) = 
(0.455, 16.8), is very coincident with the result in Rcf.[5] 
, indicating the reliability of the method for the A ~ 80 
mass region. 

In the following, the progress that used to determine 
the deformation of a nuclei will be described detailly 
with the example of 80 Kr. From the fact that the 7- 
ray energy within a superdeformed bands in A ~ 80 is 
much higher than that in A ~ 160, we get that the su- 
perdeformed nuclei in A ~ 80 rotates much faster than 
that in A ~ 160 because the rotational frequency, oj, 
is approximately the half of the 7 -ray energy. Thus, 
when we predict the shape of a nucleus in A ~ 80, 
the oj must be larger. In this paper, we fixed the oj 
as O.lo^o, where the ojq = &l/yA MeV. In the three- 
dimensional calculation, the £4 from —0.04 to 0.10 is di- 
vided to 11 points. The total routhian energy in each 




e 2 cos(y+30) 



FIG. 1: The shape determination of 80 Kr. (a) the counter 
plot of total routhian. The unit of the number in the surface 
is MeV. "+" indicates a local minimum which deformation 
is (£2,7) = (0.393,28.8°). (b) the counter plot of e 4 . The 
hexadecapole deformation at symbol "+" , which has the same 
position as "+" in (a), is 0.030. 



TABLE I: Predicted TSD nuclei between 70 < A < 90 for 
rotational frequency of O.lhuio under the selection of £2 > 0.35 
and 10° < 7 < 50° 



Nuclei 


£2 7 


£4 


Nuclei 


£2 7 


£4 


72 Ni 


0.448 12.6° 


0.048 


80 Kr 


0.394 28.6° 


0.030 


74 Ni 


0.393 21.3° 


0.023 


86 Zr 


0.471 19.0° 


0.043 


76 Zn 


0.404 21.3° 


0.029 


88 Mo 


0.500 14.7° 


0.055 


76 Ge 


0.403 30.7° 


0.027 


90 Mo 


0.375 40.0° 


0.037 


78 Se 


0.396 32.7° 


0.027 


90 Ru 


0.483 23.1° 


0.045 


80 Se 


0.351 36.1° 


0.022 









(£2 cos(7 + 30°), £2 sin(7 + 30°)) point will be minimized 
with respect to the corresponding 11 points. Fig. 1(a) 
shows a contour plot of the total routhian surface in 
which each point corresponds to the same oj but different 
£4. In the Fig. 1(a), there is a local minimum marked by 
"+" which has the deformation (£2,7) = (0.393,28.8°). 
The hexadecapole parameter £4, corresponding to the 
local minimum in Fig. 1(a), is determined by Fig. 1(b) 
which is the counter plot of £4. The £4 in each grid 
point in Fig. 1(b) is got from the minimization of the total 
routhian against £4. The symbol + in Fig. 1(b) which cor- 
responds to the minimum in Fig. 1(a) has the £4 as 0.030. 
Thus, the deformation of 80 Kr at oj = O.lwo has been 
determined as (£ 2 ,7,e 4 ) = (0.393, 28.8°, 0.030). During 
the calculation, we do not add the quasi-particle energy 
(the last item in Eq.©) to the total routhian energy be- 
cause for so high rotational frequency, one or two pairs 
of particle has been broken and it has been automati- 
cally included in the rotational energy part (see Sec. IV 
for details). 

Following the step described above, we analyzed all 
of the (3 stable even-even nuclei between 70 < A < 90, 
totally 64 nuclei. The predicted TSD nuclei are listed on 
Table I. In this paper, we call the deformation with £ 2 > 
0.35 and 10° < 7 < 50° as triaxial superdeformation. So, 
in Table 1, only the TSD nuclei under the deformation 
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FIG. 2: The prediction of triaxial superdeformed nuclei. 
The solid circle, open circle and cross represents the pre- 
dicted triaxial, axial, non- superdeformed nuclei, respectively. 
It is shown that most of predicted TSD nuclei are located in 
N = 44, 46. 

condition, £2 > 0.35 and 10° < 7 < 50°, are listed. 

The location of the predicted nuclei among the nuclei 
between 70 < A < 90 is shown in Fig. 2. In this figure, 
solid circles represent the predicted TSD nuclei, open cir- 
cles axial SD nuclei, cross symbols the nuclei in which 
we did not find superdeformation. A obvious regular 
in this figure is that when N = 42, 44, 46, most of the 
nuclei have superdeformation. Especially, for nuclei of 
N = 44, 46, most of them have triaxial superdeformation. 
Apparently, the neutron properties control the formation 
of axial superdeformation and triaxial superdeformation. 
How and why do the neutron properties control the for- 
mation mechanism of TSD nuclei? 



IV. FORMATION MECHANISM OF TSD 
NUCLEI IN A ~ 80 

In Fig. 2, it is obvious that the neutron numbers of the 
most predicted TSD nuclei are 44 and 46. Only 90 Mo 
is an exception. This phenomenon indicates that the 
neutron property governs the formation of TSD nuclei. 

In order to discuss the mechanism detailly, the 80 Kr, 
predicted to have triaxial superdeformation, is selected. 
According to the Eq. © , most part of the total routhian 
energy, E Iot , S 8 heib £pair and -E S um(= E rot + -E s hcii + 
-Epair) are plotted in Fig. 3. Fig. 3(a) and Fig. 3(b) show 
the TRS elements of proton and neutron in 80 Kr, respec- 
tively. The energy scale of the contour lines is 0.28MeV. 

In Fig.3(al), the surface of proton rotational energy is 
fiat and the deformation of the local minimum is small. 
Therefore, the proton rotating energy cannot affect the 
formation of TSD shape. In Fig.3(a2), although the pro- 
ton shell correction energy has two deep minimums, the 
proton pairing correction energy, shown in Fig.3(a3), has 
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FIG. 3: The formation mechanism of 80 Kr. (a) and (b) 
are for protons and neutrons, respectively, (al), (bl) is for 
rotating energy, (a2), (b2) for shell correction energy, (a3), 
(b3) for pairing energy, (a4), (b4) is the sum of the previous 
three item. It is shown that neutron rotating energy plays a 
key role in the formation of predicted TSD 80 Kr. See text for 
details. 



two high peeks near the minimum in Fig. 3 (a) and can- 
celed the minimum of shell correction energy. The sum 
of the three type of energy, E aum r p ), shown in Fig.3(a4), 
is fiat in the center part of the surface. This means that 
the proton properties in 80 Kr is helpless to form TSD 
nuclei. 

However, the neutron properties, shown in Fig.3(b), 
is different from Fig. 3 (a). The neutron rotational en- 
ergy, shown in Fig.3(bl) decreases sharply with increas- 
ing large £2 deformation and therefore has a strong 
driving effect towards large elongation deformation. In 
Fig.3(b2), the neutron shell correction has two mini- 
mums but were canceled by the pairing energy shown in 
Fig.3(b3). Thus, the driving force in large quadrupole de- 
formation is remained. Summing the three neutron parts 
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FIG. 4: The counter plot of the sum of liquid drop energy, 
neutron rotating energy, shell correction energy and pair cor- 
rection energy. This plot shows the local minimum marked 
by "+" which is close to the local minimum in Fig. 1(a). 

of total routhian energy, E sum t n \ } we obtain Fig.3(b4). 
This figure is most similar to Fig.3(bl), having small 
driving effect to spherical deformation and large driving 
effect to large quadrupole deformation. This is very im- 
portant for 80 Kr to form TSD shape. Fig. 4 shows the sum 
of liquid drop energy and E sum / n y A large quadrupole 
and triaxial minimum appeared on this surface. Since the 
£ sum ( p ) is flat in this region, the minimum shown in Fig. 4 
exists also in the total routhian surface, see Fig.l. In the 
formation of TSD shape, the rotational energy plays a 
crucial role. Because the neutron shell correction energy 
also decreases sharply in large deformation, it also has 
the additional role to form TSD shape. 

To confirm that the reason is also effective for other 
nuclei in A ~ 80 region, we also analyzed 78 Se which 
is predicted to have TSD shape and 86 Kr which is pre- 
dicted to have no TSD shape. The results supported our 
analysis of the formation mechanism of TSD nuclei that 
the rotational energy plays a key role and neutron shell 
energy plays an additional role in the formation of TSD 
nuclei. 

It has been pointed out that the rotational energy 
is the difference between the expectation value of H u 
(Eq.lJSJ) with and without rotation. When the rota- 
tional frequency is high, the pair of protons and/or neu- 
trons will be broken and their angular momentum align- 
ment will affect the rotational energy. In order to see 
the broken pair of protons and neutrons, the calculated 
q.p. routhians are presented in Fig. 5(a) for protons and 
Fig. 5(b) for neutrons. In Fig. 5(a), [431]3/2 orbit crosses 
with the [440] 1/2 orbit at w = 0.082?iw for protons, 
while in Fig.5(b), [420] 1/2 orbit crosses with the [413]7/2 



orbit for both a = ±5 in ui = 0.09^0 for neutrons. 
Therefore, When 80 Kr rotates at to = Q.lOhujQ, there is 
one proton pair and two neutron pairs are broken and 
this will take effect, on the rotational enererv. 
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FIG. 5: The quasipartical energy for protons (a) and neutrons 
(b) in 80 Kr. The following convection is used: solid lines: 
(tt — + , a = +|), dotted lines: (tt = +, a — ~ ^), dash-dotted 
lines: (tt — — ,ct = +5) and dashed lines: (tt — — , a = —5). 

Compared with the analysis of TSD nuclei in A ~ 160 
region, the quadrupole deformation parameters of pre- 
dicted TSD nuclei in A ~ 80 region are larger than those 
in A ~ 160 region. And also, the formation mechanism 
is different between the two regions. In A ~ 160 region, 
the neutron shell correction energy controls the forma- 
tion of TSD nuclei, while in A ~ 80 region, rotating en- 
ergy, which include the effect of q.p. angular momentum 
alignments, controls the formation of TSD nuclei. 



V. SUMMARY 

In summary, fixing the rotational frequency lu as 
0.1htoo,ll nuclei are predicted to have triaxial superde- 
formation under the condition of £2 > 0.35 and 10° < 
7 < 50° by the three dimensional TRS calculation in 
A ~ 80 region. Most of these TSD nuclei are located in 
the N — 44, 46 region, only 90 Mo is a exceptional. By 
analyzing the formation mechanism of TSD nuclei in the 
A ~ 80 region, we find that the neutron rotational energy 
which includes the contribution of quasiparticle angular 
momentum alignment plays a key role to form TSD nu- 
clei and the neutron shell energy plays additional role. 
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